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We present a reliable nonperturbative calculation of the QCD correction, at leading order in the 
electromagnetic coupling, to the anomalous magnetic moment of the electron, muon, and tau leptons 
using two-flavor lattice QCD. We use multiple lattice spacings, multiple volumes, and a broad range 
of quark masses to control the continuum, infinite-volume, and chiral limits. We examine the impact 
of the commonly ignored disconnected diagrams and introduce a modification to the previously 
used method that results in a well-controlled lattice calculation. We obtain 1.513(43) • 1CP 12 , 
5.72 (16)-10 -8 , and 2.650 (54)-10 -6 for the leading-order two-flavor QCD correction to the anomalous 
magnetic moment of the electron, muon, and tau, respectively, each accurate to better than 3%. 
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INTRODUCTION 

The experimental [T| and theoretical [5] determinations 
of the anomalous magnetic moment of the muon have 
both reached an accuracy that is better than six parts 
per million. This high precision reveals a discrepancy of 
over 3 standard deviations (3<r), which raises the possi- 
bility of physics beyond the standard model. However, 
the dominant error in the theory computation is due to 
hadronic effects that are currently not calculated but are 
instead either separately measured or simply modeled. 
This obscures the significance of the 3a effect and makes 
it difficult to improve the accuracy of the Standard Model 
calculation. 

In this Letter, we present a reliable lattice QCD de- 
termination of the leading-order hadronic correction for 
the muon, a{j vp , which is the single largest source of error 
in the theory calculation of a M . Additionally, we calcu- 
late the leading-order corrections a^ vp for the electron 
and a^ vp for the tau, achieving an accuracy of better 
than 3% for each. This was accomplished by introduc- 
ing a modification of the existing method that results 
in a significantly more well-controlled calculation. After 
examining all sources of systematic error and perform- 
ing our own extraction of the two-flavor contribution to 
the experimental measurements, we find agreement for 
all three charged leptons in the standard model. 

Our current computation is performed in two-flavor 
QCD, but the technique presented in this work is readily 
generalized to a realistic four-flavor calculation that is al- 
ready under way [3] . The precision of our calculation and 
the prospects for improving it demonstrate that lattice 
QCD can realistically provide a first-principles determi- 
nation of the leading-order hadronic contributions to the 
magnetic moments of the standard model leptons. 



LEADING-ORDER HADRONIC CORRECTION 

The anomalous magnetic moment ai of a lepton I can 
be written as a perturbative expansion in the electromag- 
netic coupling a. Contributions from QCD first occur at 
the order a 2 and can be written as [I] 

f°° 1 

a hvP = a 2 / d Q*J- W (Q*/ m f)Il R (Q 2 ) , (1) 

Jo V 

where mi is the mass of the lepton, Q is the Euclidean 
momentum and w(Q 2 /mf) is a known function. The 
combination II/? (Q 2 ) = H(Q 2 ) — 11(0) is the renormalized 
hadronic vacuum polarization function n(Q 2 ), which is 
defined shortly. The weight function w(Q 2 /mj) vanishes 
as (Q 2 )~ 2 for large Q 2 . This ensures that the integral 
above is dominated by the low Q 2 region, making it clear 
that a^ vp must be evaluated nonperturbatively. 

EXPERIMENTAL DETERMINATION 

The electron and muon magnetic moments have been 
measured in dedicated experiments [TJ [5]. To compare 
to the standard model prediction, the leading-order 
hadronic correction is determined by using unitarity and 
causality to relate the expression in Eq.[T]to 

r°° 1 

a l ; vp = a 2 ds-w'(s/m 2 )R{s). (2) 
Jo s 

Here w' is another known weight function and R(s) is the 
ratio of the hadronic cross section er(e + e~ — >hadrons) to 
the leptonic cross section er(e + e~ — >/i + /i~). The determi- 
nation of R{s) relies on the results of many experiments, 
and the integral in Eq.[2] has been evaluated by several 
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groups, most recently [H [SHE]- Additionally, there are 
higher-order corrections, including the so-called light-by- 
light contribution, which is difficult to measure and is 
modeled instead. 

Our calculation is performed in QCD with only up 
and down quarks, so we need to extract the two-flavor 
contribution to aj 1 ^. Inevitably, this introduces some 
ambiguity. For the purposes of comparing to our cur- 
rent two-flavor calculation, we adopt the simple proce- 
dure of rescaling the contribution to the integral in Eq.[2] 
from the energy regions between quark thresholds by the 
value of ^2 j Q 2 , where the sum runs over only the active 
quark flavors for that region and the electric charges of 
the quarks are eQ f . This neglects the very small changes 
due to the running of the QCD coupling, it ignores small 
off-diagonal contributions proportional to QfQf, and 
it disregards any complications at the flavor thresholds. 
These are all caveats that we must accept in the current 
comparisons but that will be eliminated in our ongoing 
four-flavor computation. 

Using the results from 9 , 1U| , we extract the two- flavor 
contributions to a t vp along the lines just described, giv- 
ing <££ a = 1.547(36) • 1CT 12 , a^ 2 = 5.660(47) • 
1(T 8 , and a^P- c * 2 = 2.638 (88) • 10" 6 . The errors result 
from propagating just those of OHOj. The systematic er- 
ror due to extracting the two-flavor contribution is likely 
larger than these uncertainties and must be taken into 
consideration when comparing our calculation to these 
estimates. 

LATTICE QCD CALCULATION 

The leading-order hadronic correction, a t vp , is the or- 
der a 2 contribution in a perturbative QED expansion of 
ai but it must be treated nonperturbatively in QCD. To 
this order in the QED coupling, the QCD corrections 
only modify the photon propagator. These contributions 
can be formally summed to all orders giving the hadronic 
vacuum polarization tensor, 

rWQ) = Jd 4 Xe lQ - x (n\Tj^x)j,(o)\n) . 

The current J M = J2f Q fQfl^Qf * s the hadronic com- 
ponent of the electromagnetic current and the sum runs 
over all relevant quark flavors. The current is con- 
served, consequently this correlation function satisfies a 
Ward identity that allows us to write Tl^ in terms of a 
single scalar function of Q 2 as 

IWQ) - {Q^Qu ~ Q 2 S^)Tl(Q 2 ) ■ 

Note that both n Ml/ and II are calculated directly in Eu- 
clidean space without any analytic continuation. 

We use standard lattice QCD techniques to 
calculate the vacuum-to-vacuum matrix element 



(il\TJ^(x) J^(y)\il). The functional integral that is 
implicit in this correlation function is evaluated stochas- 
tically using the results of the European Twisted Mass 
Collaboration [IT]. We have used two lattice spacings, 
a = 0.079 and a — 0.063 fm, to examine lattice cutoff 
corrections. Two finite-volume studies were performed 
to check for finite-size effects. The up and down quark 
masses m q are equal and are parametrized in terms of the 
pseudoscalar meson mass mps, with m q oc m PS in the 
chiral limit. As is common, we use heavier-than-physical 
quark masses and then take the limit as mps approaches 
the physical pion mass m w . This was done by studying 
the dependence on mps over the range from 650 to 
290 MeV. The so-called disconnected diagrams, ignored 
in all previous calculations, were included for almost half 
of the ensembles used in this work and are accounted 
for as a systematic error along with those from the 
continuum, infinite- volume and physical quark-mass 
limits. The additional details are standard and deferred 
to a later publication. 

Apart from variations in how the lattice calculation 
of II(Q 2 ) is matched to a smooth function, the method 
used so far in all calculations [U IT2TIT4] proceeds by nu- 
merically integrating Eq.[l] directly to form a^ vp . In our 
calculation, we parametrize T1(Q 2 ) over the entire range 
of Q 2 that is determined from our lattice computation. 
The presence of the lattice cutoff and the restriction to fi- 
nite volume induce an ultraviolet cutoff Q^ v proportional 
to 1/a 2 and an infrared cutoff proportional to 1 /L 2 . Ex- 
trapolating the functional form for T1(Q 2 ) to Q 2 = [15) . 
we numerically evaluate the integral from Q 2 = to 
Q 2 = Q 2 V . This is done without any use of perturbation 
theory, giving a completely nonperturbative evaluation 
of af vp . The systematic error caused by extrapolating to 
Q 2 = is eliminated as L is taken large and the error due 
to truncating the integral at Q 2 V is removed as a goes to 
zero. Thus both effects are automatically accounted for 
as part of the corresponding systematic errors. 

Our results for the muon using this method, which we 
deem the standard method, are shown as the lowest set 
of points in Fig.[T] Consistent with all other lattice calcu- 
lations of a|7 p , we find that the values calculated at mps 
heavier than m w are significantly lower than the experi- 
mentally measured value and apparently rise rapidly only 
when mps approaches quite near the physical value m n . 

We attribute this behavior to the contributions of 
the lowest-lying vector mesons. The rho, omega, and 
phi mesons account for over 80% of the fully measured 
ajj vp [9]. Any description of the vector-meson contribu- 
tion to aj lvp will depend on the mass my and a variety of 
dimensionless couplings. Without loss of generality, we 
focus on just those models based on my and the electro- 
magnetic coupling gv with (fijJ^lU, e) = m^gye^/^/2. 
The coupling shows up quadratically and dimensional 
analysis then results in a vector-meson contribution of 
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FIG. 1: Comparison of methods for a 1 ^ vp . The upper set of 
points are the results for a£. vp using H = my, the middle 
set use H = fv, and the lower set correspond to the stan- 
dard method, formally H = 1. The two lines are linear ex- 
trapolations of a^ vp and the curve is the phenomenological 

extrapolation of ajl vp . The three methods agree at the phys- 
ical point, denoted by the dashed line, and agree with the 
estimated two-flavor contribution to the experimental value. 
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FIG. 2: Phenomenological model for my- A model function 
is used to parametrize both our lattice calculation of my and 
the PDG value of the physical m p . This model is only used 
to illustrate the difficulties in the standard method. 



9vf( m l/ rn v)- Additionally, f(mf/m v ) should 



hvp 
a l,V 

vanish for mi — ¥ and my —¥ oo. Thus on rather 
general grounds we expect a^ V y rs Cg v mf jm v with a 
model-dependent constant C. 

These expectations can be combined with our lattice 
calculation of my and gy. As shown in Fig.[2j we find 
that my decreases moderately with decreasing mps but 
the values from our calculation are still rather high com- 
pared to the experimental result m p . Thus at some point 
a rapid decrease in my must occur. In contrast, gy, not 
shown but well fit by g v = 0.29(1) - 0.09(2) m 2 PS , has 
a mild dependence on mps and extrapolates smoothly 
to the experimental value g p . When combined with the 



model expectation ajy cx g v /m v , the behavior of aj7 p 
in Fig. [I] becomes plausible. The values of a^ vp are lower 
than the experimental value and vary moderately for the 
region oimps covered in our calculation. Only at lighter 
values of mps do we expect a sharp increase in aj^ vp . 

We can make these observations more precise, at the 
expense of introducing model dependence, by consider- 
ing the tree-level form for the vector-meson contribution 
a^y as given from effective field theory [13] . This gives 
a specific result for f(mf/m v ) that we combine with our 
calculation of my and gy to construct a model-dependent 
extrapolation of the results for a^ vp . Additionally, con- 
straining my to approach m p as shown in Fig. [2] gives 
the lowest-lying curve in Fig.[l] The apparent agreement 
with the physical value for a{j vp increases the plausibility 
that our explanation is correct. However, this construc- 
tion does not provide a reliable means of extrapolating 
our results to the physical m w but instead serves to il- 
lustrate the apparently strong mps dependence in the 
standard method. 

The difficulties encountered in the standard method 
can be traced to the occurrence of two distinct scales, mi 
and ray. Apart from any model, this is relevant because 
af vp is made dimensionless at the expense of introducing 
an external scale mi that is completely unrelated to the 
scales of QCD. Based on this observation, we define the 
following class of observables: 



hvp 



dQ 2 ^w((Q 2 /mf)(H^ ys /H 2 ))U R (Q 2 ) 



(3) 

where H is any hadronic quantity, understood to be a 
function of mps, and i? p h V s is its physical value. The 
natural choice for our calculation is H = my, but any 
choice produces a new modified quantity that has the 
same physical limit as a, yp . This follows simply by con- 
struction because H(mps — > ra n ) = H p h ys . The stan- 
dard method can be formally reproduced by the choice 
H = 1, but choosing a dimensionful scale has the ad- 
ditional advantage that the explicit dependence on the 
lattice spacing is eliminated. At the same time, the renor- 
malization condition that defines the physical limit is now 
given by the dimensionless ratio mi /-ff p h ys rather than mi 
alone. 

The calculation of a^ vp using H = my and H = fy, 
the vector-meson decay constant given by fy = mygy, 
are shown in Fig.[l] All three extrapolations agree with 
each other and with the estimated two- flavor contribution 
to the experimental measurement of a^ vp . The results for 
the new method show a significantly milder dependence 
on mps- This can be understood using the model con- 
siderations earlier. Specifically for H — my, we expect 
a vector-meson contribution of 



%v 



Cg v mf/m 



which only the mild mps dependence of gy now enters. 
The demonstration that H = fy results in similar im- 
provements illustrates that any quantity sensitive to my 
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the few percent level of our current calculation. 
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FIG. 3: Calculation of au vp for all three I = e, u and r. 
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The meaning of the symbols is the same as in Fig.nl We 

show the results from our improved method using H = my. 

The results are extrapolated linearly (solid line with error 

band) and quadratically (dashed line) to the physical point 

and agree with the two-flavor contribution extracted from the 

experimental measurements. 

will likely yield a well-controlled observable cij Vp - 

Without regard to any particular model or the experi- 
mental measurements, we can examine the relative merits 
of the standard and modified methods. Using the muon 
as an example, the shift between linear and quadratic ex- 
trapolations for a^ p (using H = my) is 1.7%, which is 
only a 0.6 a effect. The same results for a^ vp are 17% and 
3.5 a, indicating the presence of noticeably more curva- 
ture in the standard approach. In this case, cubic fits are 
required and give an extrapolated value of 4.1 (1.5) • 10~ 8 , 
which agrees with the more precise value of 5.72 (16)T0 -8 
that results from extrapolating <z^ vp . The same pattern 
holds for the electron and tau; thus, the lattice calcu- 
lation itself provides direct evidence that our modified 
method has a smoother approach to the physical limit 
leading to a more accurate calculation. 

Taking the modified method with H = my as our 
definition of cij Vp , we calculate all three I — e, fi and r. 
These results are shown in Fig.[3j and the extrapolated 
values at the physical point are 

*5/=2 = 1-513 (43). 10- 12 

<V 2 = 5.72 (16). 10- 8 

a rJ f =2 = 2.650 (54) -10- 6 . 

The quoted errors are due to the stochastic integration 
only. We do not find any statistically meaningful un- 
certainties due to lattice artifacts, finite-size effects, the 
extrapolation in mp g or the exclusion of the disconnected 
diagrams. At some higher precision these effects will be 
relevant, but there is no sign that they are significant at 



CONCLUSIONS AND OUTLOOK 

We have performed the first lattice QCD calculation 
of the leading-order QCD correction to the anomalous 
magnetic moments a| lvp that included dynamical quarks, 
examined lattice artifacts, checked finite-size effects and 
studied the disconnected diagrams. We examined the 
pitfalls of the standard method for calculating al vp and 
introduced a modification that creates a dimensionless 
quantity aj Vp composed of hadronic scales only. This 

quantity has the same physical limit as a^ vp but has a 
mild approach to that limit that is now well controlled. 
This allowed us to calculate the leading-order correction 
for all three charged leptons with an accuracy better than 
3%, reproducing our estimate of the two-flavor contribu- 
tions to the experimental measurements. 

The calculation was done using two-flavor QCD, which 
is the most significant systematic error. To resolve this, 
we are currently starting a four-flavor calculation. This 
will eliminate any ambiguity regarding the extraction of 
the two-flavor experimental value. When combined with 
further anticipated improvements, the modified method 
presented here should produce a result precise enough to 
replace the experimentally estimated a^ vp with a com- 
plete first-principles QCD calculation and eliminate this 
source of ambiguity in the current 3a discrepancy in a M . 
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